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Abstract 

Using the technique developed in approximation theory, we construct 
examples of exponential families of infinitely divisible laws which can be 
viewed as e-deformations of the normal, gamma, and Poisson exponen¬ 
tial families. Replacing the differential equation of approximation theory 
by a g-differential equation, we dehne the q-exponential families, and we 
identify all q-exponential families with quadratic variance functions when 
|q| < 1. We elaborate on the case of q = 0 which is related to free con¬ 
volution of measures. We conclude by considering briefly the case q > 1, 
and other related generalizations. 
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1 Introduction 

1.1 Exponential Type Approximation Operators 

C. P. May ism introduced exponential type operators as 

(1.1) Sx{f)im) = f Wx{m,u) f{u)du, 

_ JR 
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where Wx is a generalized function satisfying the generalized differential equa¬ 
tion 


( 1 . 2 ) 


dWx 

dm 


AVhA 


u — m 
v{m) ’ 


A > 0, 


and u is a polynomial of degree at most 2. Moreover May assumed that Sx is a 
positive operator and that 


(1.3) 


/ Wx{m, u) du = 1. 

Jm 


May knew the exact form of S\ in all possible cases except when v has two non- 
real zeros. May proved that S\ and certain linear combinations of it approximate 
continuous functions in the sense that limA^oo Sx{f,m) = f{m). Later Ismail 
and May m extended the approximation theoretic study to the case when v is 
an analytic and strictly positive function on (A, B), a component of {t : v{t) > 
0}. They also identified Wx when v{m) = 1-1- the general case with two 
complex roots. 

In the above mentioned work, it was observed that 

(1.4) / Wx{m,u) udu = m, 

Jr 

follow from 03 and CSJ. Hence m and v{m)/X are the mean and variance of 
Wxim,u), respectively. 

The parameter A is important in approximation theory since as A ^ oo the 
variance tends to zero and Wx becomes a unit atomic measure concentrated at 
u = m. Ismail and May m observed that the differential equation 03 has 
at most one solution which satisfies the normalization C3 and makes Sx a 
positive operator. They used the notation 

/•m 

(1.5) q{m) = J ce{A,B), g{q{m)) = q{g{m)) = m. 

Moreover Ismail and May proved that 

(1.6) Sx{f,m)=j Cx{u)exip^-xJ d9^ f (u) du 

and the function (or generalized function) Cx{u) is computed by inverting the 
Laplace transform 


/ Wx{m,u) {u — m)^du = 

Jr 


v{m) 


exp 



= / C'A(it)exp 

Jr 



du. 


The above formula is 


giA 


de = 


(1.7) 


exp 


A 


6 


Cxiu) exp(Xuz)du, 










and is valid for Re 2 G Range of q{m),m G {A,B). (Compare (2-1)]-) The 
theory of bilateral Laplace transform is in EOl. 

Ismail m considered the case when v{m) has a simple zero at an end point, 
which without loss of generality is taken as to = 0. He used the notation 


h{z) := 


( 1 . 8 ) 


;(z) 


^ = ^ym) \= — exp 
c 


h{9)d9^ , ri{^):=m-c + J 


0h{9)d0. 


He further assumed that h{z) is analytic at z = 0 and V(0) 7 ^ 0- In his notation 
W\ is a discrete probability distribution and takes the form 


/ nm ^ \ 

(1.9) Wx{m,du) = '^(l)n{X) exp y-J —-^^d9j6n/\{du), 
where {(()„:n = 0,l,...} are generated by 

00 

(1.10) exp(A?7(0) = 4>n{\) C: 

n—0 


and Sa{du) is a unit atomic measure concentrated at a G R. Ismail also showed 
that W\ in ini is independent of the choice of c G {A, B). 


1.2 Exponential Families 

Fix a positive non-degenerate cr-finite measure /i on R with the property that 

L{9) = / expiOu)p,{du) < 00 

Jr 

for all C < 6* < H. Denote 

k{ 9) = \iiL{0). 

The exponential family generated by /i is the set of probability measures 

:= {Ps{du) = exp{9u — K{0))p,{du) : 9 G (C, D)} . 

For a concise introduction, see PI Chapter 2]. Most authors take 9 from the 
largest admissible interval; ref. m restricts 0 to a maximal open interval. 

This family can be conveniently re-parameterized by the mean. Since fj, 
is non-degenerate, «(•) is strictly convex so that k'{) is strictly increasing on 
(C, D); it is also clear that k is analytic on (C, D). Let 

(1.11) A= lim kA9), B— lim k'(0). 

e^c+ e^D~ 

Clearly, k' : {C,D) {A,B) is invertible, and to = k'{ 9) = Jg^uPe{du) G 

{A,B). So for 9 G {C,D) probability measure Pg is determined uniquely by its 
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mean m G {A,B). Let ijj be the inverse function to k', i.e. k'= m and 
%1){k'{9)) = 9 for all m G {A,B), 9 G {C,D). Then the probability measures 

(1.12) W{m,du) -®) 

provide another parametrization of IF (/i) . Since 

/ uW (m, du) = m, 

Jr 

this is parametrization by the means. The variance function V : (A, P) ^ K. is 
now defined as 

V{m)= J {u —m)^W{'m,du), 

compare dEl- Notice that V(m) = Kl'^ij^irn)). It is known that the variance 
function V together with {A, B) determines /i uniquely, see Theorem 2.11], 
pn page 67], or [2H1 Proposition 2.2]. 

1.3 Exponential Families and Exponential Operators 

The connection between exponential families and exponential operators has been 
noticed in o Section 5], see also [dSL Theorem 2]. Here we give a somewhat 
more precise version of this relation that allows for parameter A > 0 thus con¬ 
necting exponential operators with dispersion models |24|. 

Suppose that a non-degenerate tr-finite measure /i with exponential moments 
of order 9 G {C, D) generates exponential family with the variance function 
V{m), m G {A,B). For natural A = 1,2,... denote by /ta the A-dilation of 
the convolution power , i.e. fi\{U) := {^ * ^ * ■ ■ ■ * fi){XU). The natural 
exponential family generated by /ta is the family of measures 

■■= {P\,e{du) = exp(-0M - ^ixi9))^J,xidu) : 9 G (CA, PA)} , 

where ka(^) = Xk{9/X). In particular, V'a(w) which is the inverse of k\{9) is 
■*/'a(w) = Xip(m) and the new variance function is 

(1.13) I4(m) = Kli^pxim)) = 

Notice that since k'x{9) = k'{9/X), the limits in (11.1111 do not depend on A. 
Parameterized by the mean, the family is 

•^(fa) = {Wx{m, du) : m G (d, B)}. 

We now verify that these measures satisfy equation 113. 

Proposition 1.1. If a positive non-degenerate a-finite measure p with expo¬ 
nential moments of order 9 G (C, D) generates the natural exponential family 
with the variance function V{m) defined for m G {A,B), then for natural X 
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measures fj,\ generate the natural exponential family W\{m,du) such that the 
corresponding integral operators 

S\{f){m)= J f{u)Wx{m,du) 

are the exponential type operators which satisfy equation El) with A = 1, 2,... 
and v(m) = V(m) for m G (A, B). 

Proof. It is straightforward to verify that El holds with v{m) = V{m), X = 
1,2,.... Since 


Sxif){m) 


f{u) exp{'ijjx{m)u - Kx{i>x{m)))p,x{du), 


differentiating under the integral sign we get 




i{du) 


= / f{u)ilj')^{m){u-K'^{ilJx{m)))exp{'ijj{m)u-Kiijjim))) fixidu) 


As K'^^tpxi'm-)) = m and tp'xi'm) = 1/n'f^tpxi'm-)) = 1/14(w) = X/V{m), (11.211 
follows. □ 

Remark 1.2. If equation El has solution S'a(/, m) for all 0 < A < 1, and 
m G {A, B) then the exponential family generated by p, consists of infinitely 
divisible probability laws. 

Proof. To prove infinite divisibility, without loss of generality we may concen¬ 
trate on fixed Wi(mo, p) G P{p). It is well know that with the range of means 
[A,B) kept fixed, T[p) = T{Wi{mo, p)), see |2I Exercise 2.12]. 

For A = 1/fc where fc = 1, 2,..., let Wx{m,du),m G {A,B) be the solution 
of II1.2II . The variance function is V{m)/X = kV{m). Denote by v the dilation 
of measure Wx{mo,du) by k. By Hl.ldll . the exponential family has 

the same variance function V (to) as the exponential family J^(1 Fi(too, du)). By 
uniqueness of parametrization by the means, Wi(too,c?m) = v*^{du), so infinite 
divisibility follows. □ 

Pronositiou ll .fi shows that the celebrated result ini Section 4] can be derived 
as a consequence of 123 Theorem 3.3]; the latter paper contains also several 
cubic variance functions and other interesting examples. Another interesting 
result 123 Proposition 4.4] is a consequence of pn Theorem 3.8]. 
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1.4 Notation 


We shall follow the terminology in m for hypergeometric functions, namely 


that 

n—1 

(a)„ := 1, (a)„ = ll(a + j), 

(1.14) 

3=0 

‘ ‘ 1 ■= \ (■=)»"! 


The modified Bessel functions are d 


(1.15) 


OO 


Iu{z) := ^ 

n—0 


{zl2Y+'^'^ 
nl r(n + 1 / + 1)' 


K„{z) 


^ I-u{z) - I„{z) 
2 sin(7ri/) 


The Lagrange expansion theorem d (L), page 145] says that if f{z), (j){z) 
are analytic in a neighborhood of z = 0, <('(0) Y 0 ^ := m/ip{m) then 


(1.16) f{m{O) = f{0) + J2- 


n—1 




dx''‘ 


. x—0 


By l(^a,b)(,u) we denote the indicator function of (a, b). 
Occasionally, we also use the g-notation 


: = 

(^5 Q)qO ■ 

(o-l, 0-2; ■ ■ ■ 5 Q)oo ■ — 



n 


with the usual conventions [0]^ = 
from [T7| . 


11(1 

k^O 

OO 

11(1 

(ui , Q)oO (0(2 ; Q^)oO • ■ ■ (^m ; Q^)oO ; 

1 q 

W.PI.- 

Mg'- ^ (g;g)fc(g; q)n-k 
[n-k]g\[k]ql {q-,q)n 

0, [0]q! = 1. Most of this notation is taken 


2 Examples of Variance Functions 

2.1 e-Deformations of Quadratic Variance Functions 

Letac and Mora EH page 3] raise the question of classifying exponential families 
with variances functions of the form 

(2.1) P{m) + Q{m)Y R{m), 


6 















where P,Q,R are polynomials of degree at most 3,2,2 respectively. Letac 
|28l page 74] initiated the study of variance functions 12.111 when P is a mul¬ 
tiple of R. The latter class was investigated by Kokonendji m who also 
gave an excellent overview of other known cases. Kokonendji m used prob¬ 
abilistic techniques to investigate variance functions in the Seshadri’s class 
V{m) = R{m)P{yJR{m)). This section further advances the investigation 
of the variance functions 123. 

We use Proposition o to identify certain exponential families Te with the 
variance function of the form 

(2.2) K(m) = (am^-I-6m-I-c)\/r+Tr^, e > 0. 

These are e-deformations of the quadratic variance family Pq analyzed in |23| 
and m From Mora’s theorem Theorem 2.12], as e —*■ 0 while {A, B) is 
fixed, the corresponding probability laws in weakly converge to the respective 
laws in Tq. 

We also give two examples of the functions eu which are not the variance 
functions. 


2.1.1 Continuous Exponential Families 

In this section we consider the following continuous ^-deformations: 

(i) e-Gaussian family V{m) = (1 -b em^)Vl + em^, 

(ii) e-gamma family V{m) = m^Vl + em^. 

The first case 0 gives an infinitely divisible family introduced in |S], see EZl 
Example 2.5] and |241 Exercise 3.2]. 

Theorem 2.1 (Kokonendji [27] j. For A > 0, e > 0, the exponential family 
with the variance function 

V{m) = y( 1 + -b enp, m € K. 

A 


consists of the infinitely divisible probability laws with the densities 


(2.3) exp I — 


A / 1 -b ume 


e VvTTenP 


- 1 


A 


W 


■KEpl + eu^ e 


Before we give a proof of Theorem 12.1 1 we show how we give a formal argu¬ 
ment. In the present case we have 


c = 0, q(rn) = 


m 


vT 


■ em 


9{z) = 


a/I — 


9d9 

v{9) 


= 1/^- 


:\/Y 


■ em^ 


Now after z i-^- ^zl\ u i-^- ufy/e, and A i—> Ae becomes 


(2.4) exp — = / exp{uz) C\{u)du/^/e, 

'a ^ Jr 


Re z G (—A, A). 
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If we know that the left-hand side of the above equation is a bilateral Laplace 
transform we can use the inversion theorem, Theorem 5a on page 241 of Widder 
12 a §6.5], and see that 


2 plOO 

e~^C\i—u)/y/e=- — / exp(—exp(uu) dv 

2^ ro° _ 

= — exp (—\/+ v'^) cos{uv)dv. 

^ J —00 

Formula (26), page 16 of ^21 implies 


(2.5) 


Ca(u) = 


■Vl + U 


Proof of Theorem 2.2. We verify directly. With the above C\{u) the right- 
hand side of (ITl) is 


( 2 . 6 ) 


Ae^ /■°° _ , Xi(Av^rT^) 


cosh('U2;) ■ 


•\/l + vf 


du. 


In view of (7.2.40)] 


we apply (7.14.46)] and conclude that the expression in (12.611 equals the 
left-hand side of (ITll . Substituting back the original values of A, z, u we get 

(12. at . □ 


We now consider case ®, which yields the infinitely divisible distributions 
from Hi page 58]. 

Theorem 2.2 (Letac {281 page 46, Example 8.2]). For X > 0, s > 0, the 

natural exponential family with the variance function 

1 /( 771 ) = —— vT+ern^ 

A 

defined on m > 0, consists of the absolutely continuous infinitely divisible prob¬ 
ability laws 


(2.7) 


1 -I- + cm? 

^/em 


A ^ / Au(Vl + em2) \ 

- Ix{y/eu)exp -l(o,oo)(M)du. 


Proof. We choose c = 1/x/e and apply 



e 

v{e) 


de = In 


/ mx/e{l -I- V2) \ 

y 1 -I- Vl + £777^ j 





VTTern? 

m 
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Therefore gives 


= r Cx{u)e^^^'^eicp (-Xu{\/1+ em^)/m\ du. 

y 1 + V1 + em^ J Jo ^ ^ 

To invert the above Laplace we set w = {\/l + emJ^jm so that m = Xj'JvP- — e. 
Thus for w > ■’/e we need to invert 

Cx{u) exp(-Auu;) du = + V2)^ (w + \/w^ - e). 

We use (28), page 240 in m to invert the above Laplace transform and establish 

(ITTIl . □ 



2.1.2 Discrete Exponential Families 

In this section we consider the following cases: 

(i) the £-deformation of the Poisson family V (m) = m^/l + em?, 

(ii) the discrete £-deformation of the Gaussian family V(m) = Vl + em^. 

We first consider case (jij). In this case -B = +oo and we choose c = xfs- It 
is a calculus exercise to derive 


J hie) de = In 
Okie) de = In 


1 + x/2 \ 

1 + vTTn? j 

m + vTT m? 


1 +V2 


+ 1 — m. 


Hence 


(2.8) ?(m) = 


{l + \/2)x/em ,, 1 , f Jem + Vl + em? 

i + vi + .m^ ’ ’'(«’")) = c?-TTT!— 


With C(w) = ^(to)/( 1 + v^) it follows that m = so that 

(2.9) A7?(^(m)) = In - ln(l + V2)^/^. 


A simple calculation shows that 

A oo 


i + C 
i-C 


= E 


(A). 


2 


n—0 


—n, A 
—A — n + 1 


-1 r- 


This proves the following theorem. 
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Theorem 2.3 (Letac |28l pg 98, (3)]). For A > 0, e > 0, the exponential 
family with the variance function 

V(m) = — \/l + 

A 

defined on m > 0, consists of infinitely divisible discrete probability measures 


(^■s/em + vT" 
( 2 . 10 ) 


■ em^ 


-A /^/s 


E 

n—0 


(A /^)r 


2 F 1 


—n, A/-y/e 
—A/— n + 1 


y/em 


1 + a/1 + 


5n/\{du). 



We now consider Case dni). This is again a known case: 121 Exercise 3.15] 
gives an answer in terms of the compound Poisson law, 123 Example 2.6] writes 
the answer in terms of dk(,^/£)5k{du). We remark that this is an example 

of a discrete indehnitely divisible natural family to which [221 Proposition 4.4] 
or (201 Theorem 3.3] cannot be applied. 

Theorem 2.4 (Letac [281 page 100, (8)]). For A > 0, e > 0, the natural 
exponential family with the variance function 


V{m) 


-a/T" 


■ £m^ 


defined on m > 0, consists of infinitely divisible discrete probability measures 


( 2 . 11 ) 


„--|vT+em^ 


E 7^^ E (l) d^2k-n)V^/x{du). 

'n=n ' ■ L—n ^ ^ 


n—0 ^ ’ k—0 

Proof. We choose c = 0 and apply 

dO ln(y/£m + a/1 + £m^) 




v{e) 

g{z) = q-\z) = 


sinh(2;y/£) 


v{e) 


dO = 


a/1 + — 1 


Therefore gives 

/ exp(AM 2 ;)C'A((iu) = exp (A(cosh(A/£z) — l)/£) = ^ g-A/e ^ cosh (^z) ^ 

J - r\ ^ 


n—0 


Thus 


Cxidu) =e {l)h 2 k-n)VI/xidu) 


„ (2£)"n! 

n=0 ' ' k=0 

is just the compound ^ (5 /^ + (5 ^/a)"P oisson law. Using the transform 

equation (11.611 we establish 112 .1111 . □ 
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2.2 A Rational Variance Fnnction 

Letac and Mora page 15] indicate that for pj > 0 the variance function 


{I - m/pi){l - m/p 2 ) ■ ■ .{I - m/pk) 

corresponds to a discrete infinitely divisible exponential family which is difficult 
to determine explicitly. Here we consider v{m) = to/(1 — m) which by dilation 
answers the question for fc = 1. 

In this case ^ and p of la with c= \I2 are 


^(m) = 2\/eme exp {p{^{m))) = exp = exp 2 ~ g 

With (l){z) = e*/ (2^), f{m) = exp (—^{m — 1)^ + A/8) in 11.1611 we conclude 
that 


( 2 . 12 ) ^ 


n—1 


2ngn/2^| 


pri-l J 

^e"'-exp(-A(x-l)V2) 


x=0 


n— 1 


= P-3A/8 , A/8 'S 

2”eV27i! 
ra=l k=0 


E 


n — 1 


^n—l—k 


dk+l 


exp [—X{x — 1 )^/ 2 ) 


- x—0 


For a > 0 we have 


dx^ 


rc=0 




where 


Hn{x) 


[n/2j 

E 


1=0 


JzE!E^»-2, 

j!(n-2j)! 


are Hermite polynomials. Therefore (trrni) gives the following. 


Theorem 2.5. For A > 0, the natural exponential family with the variance 
function 


V{m) 


m 

A(1 — m) 


defined on 0 < m < 1, is generated by the infinitely divisible discrete probability 
law pLx{du) = (t^n{X)5n{du) with 


(/o(A) := exp(-3A/8), 

(2.13) e-3V8 

We note that by ES Corollary 3.3] applied to the interval Mp = (0,1) we 
have (/n(A) > 0 for all A > 0. 
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Remark 2.6. One can also write 112. tail as 


4>nW = 


„-3A/8 


2n+lgn/2j^I ^ \ k 


■E 


^—k 


(fc-l)/2 


Hk+i{\/ A/2). 


Similar calculations for v{m) = m/{l + m) lead to ^ 4 ( 1 ) = —i/e/Gd, so this is 
not a variance function. 

2.3 Positivity of Wa(W', 

It is important to note that a given v(rn) does not necessarily determine a 
distribution regardless of the choice of A > 0. Ismail gave such example in Eni. 
In this section we elaborate on this example and on another example of the form 

Example 2.1. Let v{m) = m^Jl — to, to G (0,1). With c = 1/2 we find that 


^{m) = m 


1 -yr^ 


{V2-l)m . 
With C := (V2- 1)2 we have that 


, rj{^) = V2-2VT — m. 


,(0 = V 2-2 + ^. 

Therefore itrrnii becomes 

00 

(2.14) exp (a(V2 - 2) + 4ACC/(1 + CO) = E M^)C- 

n—0 

The information recorded so far is from (SO]- Comparing (4.1) and (10.2.17), 
page 189 in we see that 


(2.15) 


(j)o{X) = exp (^X{V2 - 2)^ , 

(/)„(A) = -4A(-C)"exp (A(y2 - 2)) EE (4A), n > 0, 


where Ln ^^(a^) is the Laguerre polynomial. Now 11.911 shows that W\ is a 
probability distribution if and only if (j)n{X) > 0 at the special value of A under 
consideration. On the other hand Fejer’s formula [551 Theorem 8.22.1] shows 
that L^^^Ea) is oscillatory at large n for any fixed positive A. Thus there 
is no A for which W\ is a probability distribution. This is an instance of the 
usefulness of having the parameter A. 

Example 2.2. Let us now consider the case 

v{m) = 1 — m^. 
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We take c = 0. Thus q{m) = arcsinm, g{z) = sin z, and = 1 — Vl — 

To determine C\{u) we need to invert 

exp(A(l — cosz)) = f C\{u) du, 

Jr 


for all z, Re z G Range of q{t), t G {—Trj2,7 r/2). Formula (46) page 55 of ^1] is 

poo 

(2.16) / Rrix(a) cos(xy)(ix = - 

Jo ^ 

For large p and fixed a, (19) page 88 in M is 




Therefore llTTT)ll gives 

(2.17) 

This implies 


(p2 - a2)l/4 

[ iG„(A)e“^dx = e 

Jr 


1 

^ Jr 


— p-^cosy 


(2.18) W\(m,u) = liFiA«(A) exp ^AV 1 — — Aw arcsinm) . 

Note that Ki\u[X) is real since Kp{x) = K-p{x) but it fails to be positive 
for any A > 0. Indeed, the second derivative ^ of the right hand side of HTTTh 
at 2 / = 0 fails to be negative as it equals Ae”"^. 


3 g-Exponential Families with 1^1 < 1 


Recall that for — 1 < q < 1 the ^-differentiation operator is 

,D„/)(,) := ZhLTM t„rx,^0. 

X — qx 

The g-analogue of the differential equation Ha is 


(3.1) 

This equivalent to 


Dq^rnW{m,u) 


w{m, u) 


u — m 
V{m) 


(3.2) 


ui(m, u) 


w{mq, u) 

1 + m(l — q){m — u)/V{m) 


When R(0) ^ 0 we can rescale m and u by a dilation to make R(0) = 1. Now 
lS21i has the solution 


(3.3) 


CXD 

w{m, u) = C{u) P 

n=0 


V{q'^m) 

V (m) -|- m(l — q){m — u) ’ 
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provided that the infinite products converge. 

For compactly supported measures, the following extends the notion of ex¬ 
ponential family from q = 1 to q G (—1, !)• 

Definition 3.1. A family of probability measures 

J-{V) = {w{m, u)fj,{du) : m S (A, B)} 

is a g-exponential family with the variance function V if 

(i) /J, is compactly supported, 

(ii) 0 S (A, B) and limt^o u) = w(t), m) = 1 for all u G supp(/x), 

(hi) y > 0 on (A, B) and (13.111 holds for all m ^ 0. 

Applying to both sides of f w(m,u)fi(du) = 1, from (13.111 we deduce 
that 

(3.4) J uw{m,u)fi{du) = m. 

This shows that family TiV) is parameterized by the mean. Applying Dq m to 
both sides of we deduce that 

(3.5) J {u — m)'^w{m,u)ii{du) = V{m). 

Thus V is the variance function for T{V)] compare (11.411 . 

We now show that quadratic variance functions determine g-exponential 
families uniquely. 

Theorem 3.2. If J-{V) is a q-exponential family with the variance function 

V (to) = 1 -1- am + bm^ 
and b > —1 + maxjg, 0} then 

1 -I- amq^ -\- bm'^q^^ 

1 -|- (a — (1 — q)u)mq’^ -|- (6 -|- 1 — q)m^q^^ 

and fi{du) is a uniquely determined probability measure with the absolutely con¬ 
tinuous part supported on the interval < u < -I- and 

no discrete part if of < 4b 

We remark that for b > 0 the above family of laws fi appears in m in 
connection to a quadratic regression problem. When g > 0, one could also allow 
b = —l/[N]q for some integer > 1 yielding a discrete measure p supported 
on At -I- 1 points, compare (I131l when 6 = — 1. 


OO 

(3.6) w{m,u) = 
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Proof of Theorem Vi.iH We rewrite (EH) as 

V{m) 


(3.7) 

Thus 


7 (m, u) = 


-w{qm, u). 


w{m, u) = u) 


V{m) — (1 — q){u — m)m 

1 + amq^ + 


k=0 


1 + (a — (1 — q)u)mq^ + (6 + 1 — q)nrPq 


2„2k 


from which 13. till follows by taking the limit as n ^ oo. 
We now recall that for \t\ small enough, 


(3.8) 


-j-n 

w{t,u) = V ----p„(w), 

n F 9 ’ 
n—0 


is the generating function of the monic Al-Salam-Chihara polynomials 


(3.9) UPn{u) =Pn+l{u) + a[n]qPn{u) + (1 + h[n - l]q)[n]qPn-l{u). 

This holds because the right hand side of J3.8|l satisfies lEH) , see 12]. Since p is 
compactly supported, we can integrate (EH) term by term for |t| small enough; 
we deduce that / pn(u)p(du) = 0 for all n > 1 . This determines probability 
measure p as the measure of orthogonality of polynomials Pn ■ 

Explicit formulas can be read out from Chapter 3], see also [7]- To use 
these results, we reparameterize 13.911 as follows. Let Pn{x) = a "p„(aa: + fS) 
with a = fi = ajifV — q). Then Pn{x) satisfy the three step recurrence 

(x - dq")pri(x) =Pn+i{x) + (1 -bq'^~^)[n]qPn-i{x) 


with a = 


/l-qpb+l — q 


6 = 


b+l-9' 


□ 


The technique we used in the proof of will not work beyond polynomials of 
degree at most 2. Al-Salam and Chihara 0 proved that the only orthogonal 
polynomials {pn{x)} with the generating function 


(3.10) 


^p„(a;)t” 

n—0 


oo 


n=0 


1 - axHitq^) 
1 — bxK{tq^) 


where A, H, K are formal power series with 


OO OO OO 

(3.11) A{t) = Y,ane, H{t) = Y,Kr, K{t) = Y,Kr 

n—0 n—1 n—1 

with ao/iifci ^ 0 and |a| + \b\ ^ 0 are the Al-Salam-Chihara polynomials if 
ab = 0 and the g-Pollaczek polynomials if a5 ^ 0. Theorem 13.21 corresponds to 
a = 0. The g-Pollaczek polynomials are in CH. 

(A related result appears also in |3 Theorem 23].) 
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4 Free Exponential Families 

The case q = 0 can be analyzed more directly. Since it is related to free convo¬ 
lution of measures, it is of interest to elaborate explicitly on the details. 


4.1 Free Convolution of measures 

We recall the analytic definition of the free convolution of compactly supported 
probability measures due to Voiculescu m. see also pT7l Section 2.4], ^1 Chap¬ 
ter 3]. The Cauchy-Stieltjes transform 


(4.1) 


Gfj,{z) := [ — —^i{du) 
J Z-U 


of a probability measure fj, is analytic in Iftz >0. It is known that its inverse 
G~^{z) exists for \z\ large enough. The i?-transform of ^ defined as R^{z) = 
G~^{z) — i/z plays the role of the cumulant generating function. A probability 
measure ^ is the free addictive convolution of probability measures if 

Rfj,{z) = Rf,i{z) +R^^{z). 

We write /r = /ri ffl ^2 ■ 

The free cumulants of fj, are the coefficients of the expansions 


(4.2) 


00 


= 2^ kni^J-)z^‘ 


4.2 Exponential Families with g = 0 

As previously, we consider A < 0 < B and assume that F > 0 on {A, B). 

Definition 4.1. A free exponential family with the variance function V(m) > 0 
in a neighborhood of 0 is a family of probability measures of the form 


(4.3) 


HV) := 


V(m) 


qi{du) : m G {A, B) > , 


V (jn) + m(rn — u) 
where /r is a compactly supported probability measure. 

It is easy to verify that defines a family of measures which fulfills all 
the requirements of Definition 13.11 including equation (EU with g = 0. It is 
also clear that the interval {A, B) must be chosen so that the integral (14.311 
converges. 

For the purpose of determining measure /r alone, the role of the interval 
(A, B) is insignihcant. Namely, if F is a real analytic function at 0, then ^ 
is determined uniquely by F. Indeed, since F(0) ^ 0, the Cauchy-Stieltjes 
transform EID is well defined for all real z = m + large enough, i.e. for 
all m close enough to 0, and is given by 


(4.4) 


G^{z) = 


V (m) 


16 






This determines (z) uniquely as an analytic function outside of the support 
of 

In particular, with V{m) = 1 + am + bm?, equation z = m can be 

solved for m, giving 


m = 


z — a — 


\J{a-zf -A (1 + 6 ) 

2 (1 + 6 ) 


and 

a-\- z + 2b z — \l (a — z)'^ — 4 (1 + 6) 

-- 2 ( 1 ++- -■ 

This Cauchy-Stieltjes transform appears in m (2)] in a non-commutative 
quadratic regression problem. It also appears in ^ Theorem 4], |31], and cni 
Theorem 4.3]. The corresponding laws are the free-Meixner laws 


(4.5) fj,{du) — 2'K{bv? + au + 1) ^(o-2+T+b.o+2+T+b)'^^ + + P2<JtX2 ■ 

The discrete part of /r is absent except for the following cases: 


(i) if 6 = 0, > 1, then pi = 1 — 1/a^, ui = —1/a, p2 = 0. 

(ii) if 6 > 0 and > 46, then pi = max|o, 1 — p2 = 0, and 

Ml = ± with the sign opposite to the sign of a. 

(iii) if — 1 < 6 < 0 then there are two atoms at 


Ml,2 


-a ± y/a? - 46 
26 


Pi,2 = 1 + 


- 46 + a 
26Va2 -46’’ 


This proves the free version of [221 Theorem 3.3], see also [^ Section 4]. 
Theorem 4.2. The free exponential family with the variance function 

V (to) = 1 + am + bm^ 

and 6 > — 1 consists of probability measures llOl with p given by (lOll . 

We remark that if 6 > 0 then p is infinitely divisible with respect to free 
convolution. In particular, up to dilation and convolution with degenerate 5m 
(i.e. up to ’’type”) measure p is 

(i) the Wigner’s semicircle (free Gaussian) law if a = 6 = 0; see [^ Section 
2.5]; 

(ii) the Marchenko-Pastur (free Poisson) type law if 6 = 0 and a + 0; see 123 
Section 2.7]; 
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(iii) the free Pascal (negative binomial) type law if 6 > 0 and a? > 46; see |d4l 
Example 3.6]; 

(iv) the free Gamma type law if 6 > 0 and = 46; see P Proposition 3.6]; 

(v) the free analog of hyperbolic type law if 6 > 0 and < 46; see P Theorem 

4]; 

(vi) the free binomial type law if — 1 < 6 < 0; see Example 3.4]. 

We conclude this section with the free version of Proposition o Recall 
that the A-fold free convolution is well defined for the continuous range of 
values A > 1, see m Let fi\ be the dilation by A > 1 of the free convolution 

m \ 

power 

Proposition 4.3. If a compactly supported probability measure p, generates the 
free exponential family with the real-analytic variance function V > 0 on 

{A,B), then for all X>1 measures p\ generate the free exponential family with 
the variance function V{m)/\. 

Moreover, if for every 0 < A < 1 there is a p\ which generates the free expo¬ 
nential family with the variance function V(jn)/X, then p is infinitely divisible 
with respect to the free convolution. 

Proof. From 14.411 we determine 

(4.6) R, = m, 

which determines the R-transform uniquely in a neighborhood of 0. Repeat¬ 
ing the same calculation with the i?-transform of measure p\, we get 

( m \ 

Thus 

f /a ^ 1 

J V {m) -\- Xm{m — u) 

for all jmj small enough, and p\ generates the corresponding free family US. 
The second part follows from the relation 


R 


Ux 




which proves that p is infinitely divisible with respect to the free convolution. □ 


Remark 4.4. Combining with OHll we see that the free exponential fam¬ 
ily with the analytic variance function V is defined by the unique centered 
probability measure p with free cumulants 


hn+l {p') 


■ 1 

nl dt^~^ 


V^{t) 


t=o 


n = 1 , 2 ,.... 


18 










5 ^-Exponential Families for q>l 

In this section it is convenient to set q = - with 0 < p < 1, and to use again 
auxiliary parameter A > 0. The g-analogue of the differential equation is 

(5.1) Dq^rnWx{m, du) = Wx{ni, du) 

v{m) 

As previously, from Jg^W\{m,du) = 1 we deduce by (/-differentiation that 
/jj uW\{m, du) = m and J^iu — m)‘^W\(m, du) = v{m)/X. 

With 


Ai := A(1 -p) 


we find that EU is 

(5.2) W\{m,du) = W\{pm,du) 


1 + 


Ai(u — pm)m 
v(jpm) 


We first consider (lOl when q = oo. The general case of n(0) ^ 0 reduces 
to the case n(0) = 1. Substituting p = 0 into (15.211 we get 


W\{m, du) = (1 -I- Xmu) C\{du), 

which is an analog of equation corresponding to q = oo. From this, 

it is clear that any probability measure C\{du) such that / uC\{du) = 0, 
f u^Cx{du) = 1/A determines its own (/-exponential family as long as 1 + Xum > 
0 on the support of C\{du). Moreover, it is easy to see that the only choice of 
v{m) is a quadratic polynomial v{m) = 1-1- bm — Xm? where b = la^ f u^C\{du). 
This show that g-exponential families for q = oo are not determined uniquely 
by their variance functions. 

It is plausible that non-uniqueness persists for all g > 1. For example, the g- 
Hermite polynomials {/i„(x|g)} in |22| correspond to probability measures which 
are not determined uniquely by moments. The TV-extremal solutions of the 
moment problem, pQ, are given by a one-parameter family {p{du; a) : a € (p, 1)} 
which is completely characterized by Ismail and Masson in EH, see also Chapter 
21 in |21| . Unfortunately, the construction of the corresponding exponential 
family via equation EH) led us to the family of measures 


OO 

(5.3) W\(m, du) = (1 + Ximu/q^ — Xirri^/q^^)p{du] a) 

k=0 

with negative densities. 

The non-uniqueness within the class of quadratic variance functions v(m) is 
confirmed by the following two examples. 


Example 5.1. Consider the absolutely continuous family with support in 
(0, oo) with the density 


(5.4) wx(m,u) 


(p A _ (p;p)^ sin(7rA) _ u^ ^ _ 

irm^ {p^-^;p)oo {-u{p-^ - l)/m-p)oo 
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This is the case of p-Laguerre polynomials ESI §3.21]. With 9 = 1/p, a calcula¬ 
tion verifies that 


(5.5) 


Dq^rnW\{m, u) 


TO ^(1 — p ) 


w\{m, u)(u 


Now (3.21.2), page 108 of EH shows that 


m). 


(5.6) 


w\(m, u) du = 1, 


w\{m, u)udu = m. 


(The latter integral follows also from the former by ^-differentiation and 15.511 .1 
Thus 


^ = {w\(m,u)\u>odu : m > 0} 


is parameterized by the mean. Applying again, we get the variance func¬ 

tion 


nOO 2 

(5.7) V(m)= / {u — m)'^ W\{m,du) = 

Jo \ 

with \q = . This is a continuous g-analogue of the gamma family with 

v{m) = m^. 

Example 5.2. For m > 0, consider the family of discrete measures 


W\{m, du) = wx{m, u)^{du) 
(-c,-p/c;p) 


with the density 
(5,8) 

with respect to discrete measure 

{p^;q) 


, c=(p ^ - 1)1 


fj,{du) = 


{P',P)c 


Spr^idu). 


n— — oo 


This is again related to p-Laguerre polynomials ESI §3.21]. With q 
calculation verihes that 15.511 holds. 

Now (3.21.3), page 108 of ESI shows that 


1 /p, a 


(5.9) 


/ w\{m,u) fi{du) = 1. 

Js. 


As previously, applying Dq „i to both sides of H5.9II and using 15.511 we get 


(5.10) / uwx{m,u) fi{du) = m. 

Js. 

Applying Dq^m to both sides of (OIHl and using 15.511 again, we get V{m) = 
mJjXq, compare 15.711 . Thus {Wx{m,du) : m > 0} is a discrete g-analogue 
of the gamma exponential family; it shares the variance function and the q- 
differential equation with the continuous q-analogue of the gamma exponential 
family from the previous example. 
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6 Shifted ^-Exponential Families 


The special role played by 0 in Definition Id. II is due to the fact that q-derivative 
is dilation invariant but not translation invariant. More generally, we 
consider the L-operator introduced by Hahn nsi This is a g-differentiation 
operator centered at 0 € K, which we can write as 


( 6 . 1 ) 


{eDq^^f){x) 


f{x)- f[qx +{I-q)e) 
{l-q){x-e) 


x^e, q^l. 


The usual g-derivative Dq^^ corresponds to 6* = 0. For 6* ^ 0 a dilation reduces 
all such operators to 0 = 1, in which case we use shorter notation 

( 6 . 2 ) Dq,x '■= iDq^x- 

Note that ^ ^ 

^q,x^ O5 I^q^x^ 1- 

With A < 1 < B and > 0 on {A,B), the shifted g-exponential family with 
variance function V is the family of probability measures 


J-e = {w{m, u) fi{du) : m G {A, B)} 


such that 


Dq,mw{m,u) = w{m,u) 


u — m 
V{m) ' 


In this discussion we are less restrictive than in Definition ld.il in the admissible 
range of values of 0 we include the end-points of (H, H), and we allow non¬ 
compact support for /r. Such a generalization is beyond the scope of this paper, 
so we give only one explicit example for 0<g< l,i3 = l and one for g > 1, 
H = 1. 


Example 6.1. Consider the case of the Wall polynomials, see |25l §d.20]. In 
this case we have a family of discrete probability measures 

g" 

W{m, du) = - :—Sqn(du), 

n=0 


where the density is 

(6.3) w;(to, u) = a*''“/^"‘^(ag;g)oo, a={l — m)/q. 


From (d.20.d) in we see that Jg^pi{u;a\q)w{m,u)du = 0, which implies 
/jj(l — ap — u)w(m, u)du = 0. Hence the family 

T = {W{m, du) : 0 < m < 1} 


is again parameterized by the mean, 


iW{m, du) = m. 
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From EHl (3 .20.3)] we calculate the variance function 

V(m) = m(l — m)(l — q). 

Now 

(6.4) Dq mw{m, u) = -r— -- w(m, u). 

(1 — q)m(l — m) 

Thus (16.311 defines a shifted analog of the g-Binomial family. Note that although 
equation (lOll makes sense also for g = 0, it then gives a degenerated law (5i, 
not the translation of a free binomial law. 

Example 6.2. For 0 < p < 1, let g = 1/p and consider the Al-Salam-Carlitz 
polynomials {Vn'^\x;p)}, [23 §3.25]. Let 

(6.5) W{m,du) = w{m,u)p{du) = a“ *''^(ag/u;p)ooP(dM), 

where fi{du) = /{p]p)nSp-n{du) and m = a + 1. Now with g = 1/p we 

have 

7/7 ( 7 T ) ?/ I 

Dq^rnW{m, u) = - -—2-^{1 - u(l - (ui - l)/u)}, 

(m- 1)(1 - 1/p) 

which simplifies to 

(6.6) D„ rnw{m, u) = ^ w{m, u)- — 

1 — p m — 1 

Since 123 formula (3.25.2)] implies ^^W{m,du) = 1, therefore applying 
Dq^m and taking (16.611 into account we deduce /jq uW (m, du) = m. Similarly 
V{m) = (1 — p){m — l)/p. Thus (16.511 defines the family of measures 

T = {W{m,du) : to > 1}, 

which is a shifted g-analogue of the Poisson exponential family with g = 1/p > 1. 
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